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ABSTRACT 

(.Distribution  Limitation  Statement  No.  2) 

A  hydrostatically  stable  atmospheric  model  is  necessary  to  perform  theoretical 
calculations  of  hydrodynamic  motion  in  the  atmosphere,  on  a  digital  computer. 
This  report  presents  thr  le  such  models  developed  at  the  Air  Force  Weapons 
Laboratory  for  use  in  its  hydrodynamic  computer  codes.  One  is  for  the  annu.-l 
mean  temperate  atmosphere  (45°  N  latitude);  one  for  the  annual  mean  tropical 
atmosphere  (15°  N  latitude);  and  one  for  the  summer  subarctic  atmosphere  (60° 

N  latitude).  The  models  are  presented  herein  in  tabular  form  and  as  FORTRAN 
subroutines  which  could  be  placed  directly  into  any  hydrodynamic  computer  code. 
For  a  given  altitude  (cm),  the  subroutines  return  a  pressure  (dynes/cm2), 
density  (gms/cm3),  specific  internal  energy  (ergs/gm),  temperature  (°K),  and 
(y  -  1).  The  pressures  and  densities  agree  with  tabulated  values  to  at  least 
1  part  in  108  and  temperatures  to  at  least  1  part  in  100.  The  atmospheres 
experience  an  acceleration  of  no  more  than  3  parts  in  10-  in  a  first  order 
finite  difference  scheme  with  a  zone  size  of  1  kilometer,  the  worst  case. 
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SECTION  I 
INTRODUCTION 


The  Air  Force  Weapons  Laboratory  is  engaged  in  research  or  the  phenomenology 
of  atmospheric  nuclear  weapons  detonations.  This  research  consists  of  numerical 
experiments — theoretical  calculations  of  various  yield-altitude  detonations  — 
performed  on  large  memory',  high-speed  computers.  These  calculations  are  made 
by  a  series  of  large,  complex,  computer  codes  that  enable  prediction,  from 
essentially  first  principles,  of  the  phenomenology  from  microseconds  to  minutes 
after  a  detonation. 

Three  main  codes  are  used:  SPUTTER,  a  one-dimensional  radiation  transport 
Lagrangian  hydrodynamic  code;  SAP,  a  one -uimensional  Lagrangian  hydrodynamic 
code;  and  SHELL,  a  two-dimensional  Eulerian  hydrodynamic  code  that  also  has 
the  capability  of  doing  radiation  diffusion.  The  first,  SPUTTER,  takes  the 
radiative  output  directly  from  weapons  design  calculations,  deposits  the  energy- 
in  air,  and  calculates  the  ra’iative  and  hydrodynamic  growth  of  the  fireball  to 
a  time  of  about  1  second.  It  uses  a  multifrequency  transport  scheme  involving 
typically  20  frequency  groups.  The  output  of  this  code  is  used  as  input  to 
SAP,  which  calculates  the  hydrodynamic  expansion  of  the  shock  wave  at  various 
angles  to  the  horizontal.  The  output  of  SPUTTER  is  also  used  as  input  to 
SHELL,  which  calculates  fireball  growth  and  rise  to  late  times. 

These  codes,  or  any  hydrodynamic  code  in  general,  wh’ch  perform  hydrodynamic 
calculations  in  a  real  atmosphere,  need  an  atmospheric  model  that  is  stable  in 
a  R~‘  gravity  field;  otherwise,  the  minvr  cal  integration  of  the  hydrodynamic 
equations  will  result  in  fictional  vertical  velocity  component s s 

Tins  report  describes  three  stable  atmospheric  models  developed  at  too 
Air  Force  Weapons  Laboratory.  The  first  ewdel  is  base-’  on  data  in  ref*  pence  1 
for  an  annual  mean,  temperate  atmosphere  (45*  N  latitude)  for  ul'itjdes  f ror 
se3  ieve1  to  700  kilometers.  The  secoiiu  is  based  on  data  in  reference  for  ar 
annual  mean,  tropical  atmosphere  (IS'  N  latitude)  for  altitudes  from  sea  level 
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to  90  kilometers  and  reference  1  for  altitudes  from  90  to  700  kilometers. 

The  third  is  based  on  data  in  reference  3  for  a  summer  subarctic  atmosphere 
(60®  N  latitude)  for  altitudes  from  sea  level  to  90  kilometer  and  reference  1 
for  altitudes  from  90  to  700  kilometers 

The  Appendix  contains  a  FORTRAN  subroutine  for  the  models. 
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SECTION  T I 

DERIVATION  OF  ThE  MODELS 


1 .  Conditions  and  Assumptions 

The  molecular  scale  temperature  is  the  defining  property  of  the  US  Standard 
Atmosphere,  1962  (Ref-  1).  It  is  defined  by  the  relationship 


T  (z) 
ra 


M(z) 


T(z) 


(1) 


where 


T(z ) 
M(z) 


M 

o 


molecular  scale  temperature  at  altitude  z 
absolute  temperature  at  altitude  z 
molecular  weight  of  air  rt  altitude  z 
molecular  weight  of  air  at  sea  level 


The  variation  of  the  molecular  scaie  temperature  is  defined  3S  a  series  of 
connected  s-gments  linear  in  geopetentia!  altiude  to  90  kilometers  and  in  geo¬ 
metric  a’.tltude  above  9G  kilometers  (Ref.  1).  Therefore,  we  diviat  rne  temperate, 
tropical,  and  subarctic  models  into  21,  22,  and  23  altitude  groups  respectively 
from  sea  level  to  700  kilometers  and  assume  a  linear  molecular  scale  temperature 
variation  with  altitude  within  each  altitude  group.  This  v.i  r :  it  i  on  has  tl  c 
form 


•  v(  i}  +  (Li)hi) 


where 


z 


i 


the  base  altitude  of  the  i-th  aitit-rii  group 


the  molecular  sc.  e  temperature  gradient  over 
group 
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We  3lso  assume  the  atmosphere  behaves  like  an  ideal  gas  in  which  case 


p(z) 


M(z)p(z) 

RT(z) 


(3) 


where 


p(z)  t  the  density  at  altitude  z 
p(z)  =  the  pressure  at  altitude  z 
R  ~  the  urive- sal  gas  constant 


Finally,  we  requi  '  che  atmospheric  models  to  satisfy  the  condition  of 
stability  given  by  the  Hydrostatic  ;uation 


~p(z)g(z) 


where 


g(z)  =  the  acceleration  due  to  gravity  at  altitude  z 

V" 

■  (a+z) 2 

gQ  =  the  acceleration  due  to  gravity  at  sea  level 
a  =  Che  radius  of  the  earth 


(4) 


(5) 


By  combining  equations  1  through  4,  the  hydrostatic  equation  becomes 


p(z) 


g  M  a 
o  o 


( a+z ) - 


dz 

[  V‘)“  (l‘)(M‘)] 


(6) 


Integration  of  equation  6  from  z^  to  z  results  in  the  final  form,  which  is 
Co  calculate  pressures  in  each  altitude  group. 


sIlL. 


g  H 
o  o 


(*o  ■  ■  r^)- 


-  ,(li)K)]K) 


(a+z)/ a+z 
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a+z.  T 
_ i  m 

y 

+  ( 

(a+z) 

T 

n 

w 

(7) 
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2.  Calculations 

To  calculate  the  molecular  scale  temperatures,  T^z^j,  and  the  molecular 
scale  temperature  gradients,  ,  which  define  each  atmospheric  model,  we  take 
the  pressures  listed  in  references  1,  2,  and  3  corresponding  to  the  base 
altitude  of  each  altitude  group  and  the  molecular  scale  temperatures  as  defined 
by  equation  2  tc  be  constants-  With  p^+^j,  i  a»  8Q»  R»  and 

known,  equation  7  becomes 


g  M  a- 


£n 


’(2H-l)  _ _ _ 

PW  \Wr 


2R 


[rm(zl)  -  (Ll) (a*2l)l  (z 1+1 ~ Z i ) 

pl+l)(-i) 


+  L^tn 


**i  T»(h 

l)  -  (Li/(2'+l'zi) 

*+2l+! 

N 

E 

H 

| 

(8) 


which  is  solved  by  an  appropriate  iterative  technique  to  find  that  molecular 
scale  temperature  gradient,  L^,  which  will  satisfy  equation  8  and  thus,  satisfy 
the  condition  of  hydrostatic  stability  for  the  i-th  group.  Then  the  molecular 
scale  temperature  corresponding  to  the  base  altitude  of  the  next  altitude 
group  becomes 


Tm(2i+l)  “  ^m(Zi)  +  (Li) (2i+l"Zi) 


(9) 


where  is  the  gradient  just  calculated.  Then  the  calculation  proceeds  to  the 
next  altitude  group  where  the  exact  procedure  is  repeated.  This  procedure  is 
repeated  until  the  molecular  scale  temperature  and  gradient  have  been  calculated 
for  all  altitude  groups. 

The  calculation  begins  at  sea  level  for  which  we  take  the  molecular  scale 
temperature  as  listed  in  references  1,  2,  and  3.  (This  temperature  can  be 
thought  of  as  an  "initir.l  condition."  All  other  temperatures  are  defined  by 
the  gradients  obtained  from  the  solution  of  equation  8  and  equation  9.)  Solu¬ 
tion  of  equation  8  provides  L,  which  is  used  in  equation  9  to  obtain  T2.  Using 
T2,  equation  8  gives  L2,  which  in  turn  by  equation  9  gives  T3,  and  so  forth. 
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We  repeat  the  above  iteration  foi  each  model  using  the  apnropriate  data  for 
P(2ij*  an<*  ®o  38  ^ven  ^  the  references. 

3.  Transition  Region  at  90  km  for  the  Tropical  Atmospheric  Model 

We  used  the  procedure  outlined  above  to  calculate  the  defining  properties 
of  all  three  atmospheric  models.  However,  in  the  case  of  the  tropical  model, 
the  arithmetic  error  between  calculated  and  tabulated  temperatures  increased 
in  magnitude  and  alternated  in  sign  with  increasing  altitudes  above  90  Rm. 

This  result  was  caused  by  the  poor  match  of  data  at  90  km  provided  by  refer¬ 
ences  1  and  2.  Table  I  shows  the  difference  in  data  provided  by  the  two 
references  at  90  km. 

Table  I 

DIFFERENCE  IN  ATMOSPHERIC  DATA  FOR  THE  ALTITUDE 
90  km  GIVEN  BY  REFERENCES  1  AND  2 

p(dynes/cm2)  p(gm/cm3) 

Ref  1  1.6438  3.1700  x  10“9 

Ref  2  1.8620  3.5224  x  10'9 


T  (°K) 

IP 

180.65 

184.15 


Therefore,  we  created  a  transition  region  between  90  and  110  km  to  effect 
a  smooth  transition  from  the  data  in  reference  2  (0  to  90  km)  to  the  data  in 
reference  1  (90  to  700  km).  We  did  this  by  replacing  the  two  altitude  groups 
from  90  to  100  km  and  from  100  to  110  km  by  two  groups  from  90  km  to  an  altitude, 
z,  to  be  determined,  and  from  z  to  110  km. 


We  determined  the  altitude,  z,  in  the  following  manner.  First,  the  basic 

procedure  outlined  in  section  II-2  was  used  to  calculate  the  molecular  scale 

temperatures  and  gradients  for  the  altitude  groups  from  0  to  90  km,  using  the 

data  in  reference  2.  Next,  we  assumed  that  the  isotherm*!  region  from  80.13 

to  90  km  actually  extends  from  80.13  km  to  the  altitude,  z;  that  is,  T  (80.13)  « 

m 

T  (90)  ■  T  (z)  or  L(80.13-*z)  «*  0.  Then  we  calculated  the  alti*  z,  by  itera- 
m  m 

tion  according  to  the  analysis  below. 

Since 


(10) 


we  have  an  equation  that  can  be  solved  by  an  iterative  method  to  give  the 
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altitude  z.  The  first  term  on  the  right  of  equation  10  becomes 


♦ 


4 


2 

SliL.  m  _  jgjgL.,  z-90 

p(90)  T  (90) R  (a+z) (a+90) 

Ul 


(ID 


The  second  term  becomes 


ln  eLUPI  «  _ _ goj?3 _ 

p(z)  rTn(90)  -  (L)(a+90)J^R 


li 


(90)  -  (L)(a+90)  (z-90) 


(a+z) (a+90) 


l]  (z-S 


U0  1,(90)  ♦  (10(2-90)' 

Lln  |^a+z  Tn(90) 


(12) 


The  values  p(90)  and  gQ  are  known  from  reference  2;  the  values  p(110),  a,  R, 

and  T  (110)  from  reference  1;  and  T  (90)  from  the  calculation.  Since  L  is  a 
m  m 

function  of  z 


T  (110)  -  T  (90) 
m  m _ 

110  -  z 


(13) 


equation  10  can  be  solved  by  Iteration  to  find  z. 

After  finding  the  altitude  z,  and  p(z)  by  equation  11,  we  resumed  the 
calculation  of  molecular  scale  temperatures  and  gradients  for  the  altitude 
groups  from  90  to  700  km,  using  the  method  outlined  in  section  II-2. 

Similar  differences  between  references  3  and  1  exist  at  90  km  for  the 
subarctic  atmosphere.  However,  it  was  possible  to  generate  a  hydrostatically 
stable  model  without  resorting  to  an  artificial  transition  region.  Research 
on  the  subarctic  model  will  continue  to  determine  whether  such  a  region  will 
indeed  improve  it.  In  any  case,  the  present  subarctic  model  is  adequate  for 
o-t  purposes. 
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SECTION  III 
APPLICATIONS 

When  fitting  an  atmosphere  into  a  mesh  of  cells  for  the  purpose  of  solving 
the  hydrodynamic  equations  numerically  on  a  computer,  it  is  necessary  to  define 
a  density  and  specific  internal  energy  as  well  as  a  pressure  to  the  atmosphere 
contained  in  each  cell.  Therefore,  for  a  given  altitude  z 


p(z) 


’(*i) 


-f 


(14) 


where 


and 


g  M  a" 
o  o 


T»(gi)  *  (li)(,-ii) 


T»k)  -  (tt)(,wi)l(,-,i) 

(a+z)  ^a+z^ 


+  L1tn 


ri  -  *  <  zi+l 


By  the  idael  gas  law,  the  density  at  altitude  z  is 

>  .  T  (o) 

c(a)  -  ~r  y  ,-(o) 

p(o)  T  (z) 


(15) 


(36) 


where  p(o),  r(o),  and  T  (o)  are  the  sea  level  values  for  pressure,  aensity,  and 
molecular  scale  temperature,  respectively. 

To  obtain  the  specific  internal  energy,  !(*>,  we  require  the  equation  of 
state  for  air  which  will  be  used  In  the  calculation  (e.g.,  t*  i  Doan-SicKel 
Equation  of  State  for  Air.  reference  4)  to  return  the  same  n  assure  as  the 
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•tnosph«ric  model.  Thle  involves  s  solution  by  iteration  which  proceeds  ai 
follows.  Guess  a  value  for  y.  Then,  by  the  ideal  gas  law. 


1(2) 


_ Pl*l  . 

(y-l)p(z) 


(17) 


where  p(z)  and  p(z)  are  the  values  of  pressure  and  density  at  altitude,  z,  as 
defined  by  the  atmospheric  model.  Then  these  values  for  I(z)  and  p(z)  are 
entered  into  the  equation  of  state  for  air  which  returns  a  value  for  ^ef fective 
•  J |.  Then,  by  the  ideal  gas  law,  the  pressure  is 

P(Z)  *(Yeffective“1)p(2)I(z)  (18) 

The  iteration  continues  until  that  coibbination  of  y  and  I(z)  is  found,  which 
will  make  the  equation  of  state  deliver  the  same  pressure  as  the  atmospheric 
model.  See  appendix. 
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SECTION  IV 
RESULTS 

The  annual  mean  temperate  atmospheric  model,  the  annual  mean  tropical 
atmospheric  model,  and  the  summer  subarctic  atmospheric  model  are  given  in 
tables  II,  III,  and  Iv,  respectively.  Tables  V,  VI,  and  VII  show  the 
difference  between  the  molecular  scale  temperatures  at  the  base  altitude  of 
each  altitude  group  as  listed  in  references  1,  2,  and  3  and  those  calculated 
for  the  atmospheric  models. 

The  value  for  the  altitude,  z,  used  in  the  transition  region  of  the 
tropical  atmosphere  is  97.84061  km;  that  is,  the  two  altitude  groups  making  up 
the  transition  region  extend  from  90  to  97.84061  km  and  from  97.84061  to  110  km. 

The  appendix  contains  the  models  in  subroutine  form. 
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Table  II 

DEFINING  PROPERTIES  OP  THE  ANNUAL  MEAN  TEMPERATE  ATMOSPHERE 
(45*  N  latitude) 


Molecular  acale  Molecular  acale 

Altitude  Preaaura  temperature  gradient  temperature 

*  x(cm)  (dynea/cm2)  L  (*K/cm)  Tu  (*K) 


0.0000 

1.01325 

X 

106 

-6.49291767 

X 

IO'5 

288.150000 

1.1019 

X 

106 

2.26320 

X 

105 

9.28049177 

X 

10”* 

216.604540 

2. 0043 

X 

106 

5.47487 

X 

io1* 

9.86254816 

X 

IO”6 

216.688473 

3.2162 

X 

106 

8.68014 

X 

103 

2.77080370 

X 

IO”5 

228.621170 

4.7350 

X 

10s 

1.10905 

X 

10 3 

-1.72246873 

X 

10"7 

270.704137 

5.2429 

X 

106 

5.90005 

X 

io2 

-1.95999298 

X 

IO-5 

270.616652 

6.1591 

X 

106 

1.82099 

X 

102 

-3.91697376 

X 

10“5 

252.659197 

7.9994 

X 

106 

1.03770 

X 

io1 

1.60823156 

X 

io*7 

180.575129 

9.0000 

X 

106 

1.64380 

X 

10° 

2.98166731 

X 

io-5 

180.736048 

1.0000 

X 

107 

3.00'  1 

X 

io’1 

5.02020153 

X 

io”5 

210.552722 

1.1000 

X 

107 

7.35440 

X 

IO'2 

9.97762308 

X 

io-5 

260.754737 

1.2000 

X 

107 

2.52170 

X 

io-2 

2.00108806 

X 

io- 

360.530968 

1.5000 

X 

107 

5.06170 

10”  3 

1.49589024 

X 

io- 

960.857386 

1.6000 

X 

107 

3.694  -0 

10"  3 

1.00407491 

X 

io— 

1110.446410 

1.7000 

X 

107 

2.79260 

IO-3 

6.97598503 

X 

10— 

1210.853900 

1.9000 

X 

10  7 

1.68520 

IO'3 

5.01601097 

X 

io'5 

1350.373600 

2.3000 

X 

10  7 

6.96040 

10““ 

3.98897144 

X 

10”5 

1551.014040 

3.0000 

X 

107 

1.88380 

10- 

3.31099390 

X 

10- 

18V).  242040 

4.0000 

X 

107 

4.03040 

IO"5 

2.58868496 

X 

10- 

2161.341430 

5.0000 

X 

107 

1.09570 

10~5 

1.71252931 

X 

10- 

2420.209930 

6.0000 

X 

107 

3.45020 

10"* 

1.09162418 

X 

10'5 

2591.462860 

7.0000 

X 

107 

1.19180 

10- 

... 

2700.625280 

a  -  C. 35670  *  10* 
g  *  9.80665  *  IO2 
I  •  a. 314*0  x  io7 
M  -  2.89644  x  10* 


cm/ etc1 

arg/mole/deg  K. 

gm/anle 
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Table  III 

DEFIX TNG  FX0FEXTIE8  OF  THE  AXXUAL  MEAN  TtOFICAL  ATMOS PHEKE 
(IS*  X  latitude) 


Altitude 

«(c*) 


Fracture 

(dyaee/ca2) 


Molecular  acale 
taaperatur*  (radiant 
L  CX/cb) 


Molecular  acala 
teaperature 

\  CK) 


0.000000 

1.013250 

X 

104 

-3.69962272 

10"5 

299.650000 

?. .  254000 

X 

105 

7.813300 

X 

10s 

-1.29419605 

10*- 

291.311050 

2.505000 

X 

10s 

7.586100 

X 

10* 

-6.75109163 

10"5 

288.062618 

1.657000 

X 

10* 

1.013700 

X 

10s 

4.11001824 

10*5 

193.108515 

2.211300 

X 

106 

4.043000 

X 

io- 

2.14811736 

10"5 

215.890346 

6.743200 

X 

106 

1.188600 

X 

103 

4.20150001 

10"  7 

270.278529 

5.149800 

X 

106 

7.175600 

X 

102 

-1.93134924 

10"5 

270.449362 

*  QA«‘>SV) 

X 

106 

2.535100 

s 

9  /v2 

iv 

-j.tj/iiSzo 

10"5 

254.701140 

8.013000 

X 

106 

1.100100 

X 

10 1 

1.98418410 

10"7 

184.313581 

9.000000 

X 

106 

1.862000 

X 

10® 

-9.16634469 

10"  7 

184.509420 

9.784061 

X 

106 

4.549471 

X 

10* 1 

6.35832871 

10"5 

183.790723 

1.100000 

X 

107 

7.354400 

X 

10"2 

9.74012527 

10"5 

261.104121 

1.200000 

X 

107 

2.521700 

X 

10‘2 

2.00803436 

io"- 

358.505374 

1.500000 

X 

10  7 

5.061700 

X 

10* 3 

1.44417115 

10"- 

960.915680 

1.600000 

X 

107 

3.694300 

X 

10*3 

1.05385251 

10'- 

1105.332790 

1.700000 

X 

107 

2.792600 

X 

10"  3 

6.68110562 

10"5 

1210. 718050 

1.900000 

X 

107 

1.685200 

X 

10"3 

5.15616802 

10"5 

1344.340160 

2.300000 

X 

107 

6.960400 

X 

10"- 

3.88307486 

10"5 

1550.586880 

1.000000 

X 

107 

1.88 38 

X 

10"- 

3.37947184 

10'5 

1822.402120 

4.000000 

X 

to7 

4.030400 

X 

10"5 

2.49828094 

10'* 

2160.349  300 

5.000000 

X 

107 

1.095700 

X 

10"J 

1 . 79900870 

io~s 

2410.177*30 

6.000000 

X 

10  7 

3.450200 

X 

10-* 

9.94365149 

10"* 

2590.078270 

7.000000 

X 

10  7 

1.191800 

X 

10"* 

.  .  . 

2689.514780 

a  -  6.35470  *  1J*  ca 
^  «•  9.78381  a  :0*  oe/eer2 
X  «  8.31440  *  107  er*/ao'*/<Je§  E. 

n  -  2.99644  x  101  **/•<>  1« 
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Table  IV 

OKrnme  raorotriK  or  xm  sown  sbbaktic  athosphskc 

(60*  H  latitude) 

Molecular  scale  Molecular  acala 

Altitude  Preaaure  f f aratura  gradient  teeperature 

*(c*)  (dyaea/c*2)  L  (*K/ca)  I,  (*K) 


0.0000 

1.010C 

10* 

-5.15226593 

X 

IQ"3 

287.150000 

4.9980 

X 

103 

5.4154 

10s 

-7.45507917 

X 

iO"3 

261.398975 

1.0003 

X 

106 

2.6758 

10s 

1.75752*72 

X 

10"* 

224.086304 

2.3054 

X 

10s 

3.7248 

1C* 

1.21402825 

X 

10*5 

226.380054 

3.2121 

X 

106 

9.8883 

103 

3.20422887 

ir 

10~3 

237.3876*9 

4.3237 

X 

1C* 

2.2624 

103 

5,54484005 

X 

10"* 

273.005857 

4.8303 

X 

10s 

1.2140 

103 

5.03741727 

X 

10"* 

275.81*873 

5.3376 

X 

i<* 

6.5545 

102 

-2.35321674 

X 

10"5 

278.37035* 

5.9476 

X 

i<* 

3.0773 

102 

-4.51273514 

X 

icr3 

264.015732 

7.9890 

X 

i<* 

1.2772 

101 

-1.48868732 

X 

10"* 

171.892757 

9.0000 

X 

10s 

1.7851 

10° 

3.31280965 

X 

10‘5 

170.38769* 

1.0000 

X 

107 

3.0075 

10’ 1 

6.61955253 

X 

10"5 

703.515791 

1.1000 

X 

107 

7.3544 

X 

10"2 

8.07805343 

X 

10~3 

269.711316 

1.2000 

X 

107 

2.5217 

I  o*2 

2.107*0515 

X 

10"* 

J50.49185C 

1.5000 

X 

107 

5.3617 

10~3 

1.067090*0 

X 

10"* 

982.713396 

A.  ,  - 

X 

107 

3. *943 

10"  ’ 

1.47081912 

X 

10"* 

1089.422*40 

1.7000 

X 

IQ7 

7  . 7926 

10“  3 

4.*9679987 

X 

IO"3 

1236.504250 

1.9000 

X 

107 

1.6852 

10“ 3 

6.37859718 

X 

ID"3 

1326.440340 

2.3000 

X 

107 

6.9604 

IP"* 

3.131*5277 

X 

10*3 

1581.58*230 

3.0000 

X 

107 

1.8838 

10"* 

3  9957*107 

X 

IP*3 

1800.785930 

4.0000 

X 

1C? 

4.CJ04 

10"3 

1.S38 56612 

X 

10"  3 

2200.3*0030 

5.0000 

X 

10 ! 

1,0957 

10"3 

2.522712*7 

- 

.0- 

2384.216*70 

6.0000 

X 

107 

3.4502 

10"* 

2.47207*16 

X 

’.o'* 

2636.487910 

’.GOOC 

X 

10  7 

1.1911 

10“* 

— 

26*1.208650 

a  •  6.15670  a  1C*  ce 
ge  «  9-81911  a  102  ce^eac7 
t  *  8.114*0  a  IC'  erg/aeU/d**  * 
M  -  1.8*64*  a  10*  ga^aole 
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Txblx  V 

eimueii  n.uia  ckixsoAva  urn  zabciated  •ctjc'tu* 

SCALE  TOCIIAniUS  U  THE  TDVOtATE  .THOCTTOIC  MODEL 


AltUodx 

(ca) 

txkxliul 

c«9*nt«m 

<**) 

Calculate 
txapcrxturxa 
'  *E) 

2 

Wfftrxocx 

0.0000 

208.15 

288  15 

0.000 

1.1014  I  10* 

216.65 

216.60 

-0.02i 

2.0063  i  10* 

216.65 

216.69 

0.018 

;.2162  x  10* 

228.65 

228.62 

-0.013 

6.7350  x  10* 

270.65 

270.70 

0.020 

5.2624  >  10* 

270.65 

270.62 

-0.012 

6.1541  s  10* 

252.65 

252.66 

0.004 

7.4446  x  10; 

180.65 

180.58 

-0.041 

9.u000  x  10* 

180.65 

180.7* 

0.048 

1.0000  x  107 

210.65 

210.55 

-*.0*6 

1.1000  x  107 

260.65 

260.75 

0.040 

1.2000  x  107 

360.65 

560.53 

-0.033 

i.5oe«  x  io7 

460.65 

460.86 

0.022 

1.6000  X  10T 

1110.65 

1110.45 

-0.018 

1.7000  X  107 

1210.65 

1210.85 

2.017 

1.4000  x  107 

1350.65 

1350.37 

-0.020 

2.3000  ;  107 

1550.65 

1551.01 

0.023 

i, jooe  x  ie~ 

18)0.65 

18)0.24 

-0.022 

1.0000  *  1C7 

2160.65 

2161.31 

0.032 

j.aooe  *  >.o7 

1420.65 

*  2*70.  M 

6.9000  x  I=97, 

2540.65 

2541.4* 

0.031 

7.0000  t  107 

2700-65 

2700.63 

-0.001 
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Table  VI 

oinrmscES  irtvra  calculated  a®  tabulated  moucula* 

SCALE  norOATUEU  IS  THE  TSOPICAL  ATMOSPHERIC  HDDLL 


Altitude 

(cm) 

Tabulated 

taaperatt.ea 

CE) 

Calculated 
teape raturea 
v*E> 

2 

Difference 

0.000000 

299.65 

299.65 

0.000 

2.25*000  x  10s 

286.15 

291.31 

1.304 

2.505000  x  10s 

286.95 

288.06 

0.388 

1.657000  x  10° 

193.15 

193.11 

-0.021 

2.211aiu  x  Ur 

215.15 

215.8? 

7.3** 

*.763200  i  lO" 

270.15 

270.28 

0.0*8 

5.198000  x  10* 

270.15 

270. *5 

0.111 

5.965200  x  JO* 

25' .15 

256.70 

0.217 

8.013000  x  10* 

184.15 

18*.  31 

0.089 

9.000000  x  IP4 

18*.  15 

18*. 51 

0.915 

9.78*061  x  107 

186.15 

183.79 

-0.191 

1.100000  x  107 

260.65 

261.10 

0.17* 

1.200000  x  107 

360.65 

353.51 

-0.595 

1.500000  x  10T 

960.65 

960.?? 

0.028 

1.60000C  x  107 

1110.65 

1105.33 

-0.478 

1.700000  x  107 

1210.65 

1210.72 

0.006 

1.900000  x  10’ 

1350. ',5 

13**.  3* 

-0.467 

7.300000  x  107 

1550  6  5 

1550.59 

-0.004 

3 .000000  X  107 

18TJ.65 

1822. *0 

-0.451 

*.000000  X  107 

2160.65 

2160.35 

-0.01* 

5.000000  x  10” 

2*20.65 

2*10.18 

-0.432 

6.000001  ,  l"7 

2590.65 

2590.08 

-0.022 

7.000000  *  io’ 

2700.65 

2689,51 

0.<*12 
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Table  VII 

DIFFERENCES  BETWEEN  CALCULATED  AND  TABULATED  MOLECULAR 
SCALE  TauPES/TURES  IN  THE  SUmfTIC  ATMOSPHERIC  MODEL 


Altitude 

(cm) 

Tabulated 

temperatures 

(*K) 

Calculated 

temperatures 

(-<) 

Z 

Difference 

0.3000 

287,15 

287  15 

0.000 

4.9980  x  105 

260.15 

261.40 

0.180 

1.0003  x  10* 

225.15 

224.09 

-0.472 

2.3054  x  106 

225.15 

226.38 

0.546 

3.2121  x  106 

238.65 

237.39 

-0.529 

4.3237  x  10f 

271.65 

273.00 

0.499 

4.8303  x  106 

277.15 

275.81 

-0.482 

5,3376  x  106 

277.15 

278.37 

0.440 

5.9476  x  106 

265.15 

264.0?. 

-0.-28 

7.9850  x  106 

171.15 

171.90 

0.434 

9.0000  x  10c 

171.15 

170.39 

-0.445 

1.0000  x  TO7 

210.65 

203.52 

-3.387 

1.1000  x  107 

260.65 

269.71 

3.476 

1.2000  x  107 

.>60.'' 

350.49 

-2.817 

1.5000  x  107 

960.65 

982.71 

2.297 

1.6000  x  10 7 

1110.65 

1089  a 

-1.911 

1.7000  x  10^ 

1210.65 

1236. 5C 

2.136 

1,9000  x  107 

1350.65 

1326.44 

-1.79? 

2.3000  x  107 

1550.65 

1591.58 

i.«:: 

3.0000  x  107 

1630.65 

1800.79 

-1.631 

4.000C  x  ;o7 

2’60.65 

2200.36 

1.327 

5.0009  x  107 

2420.65 

2384 .22 

-1.505 

6.0000  x  197 

2590.65 

2636.49 

1.769 

7.0000  x  107 

2700.65 

2661.21 

-1.460 
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SECTION  V 


CONCLUSIONS 


The  atmospheric  models  are  darived  such  that  the  hydrostatic  equation  is 


satisfied. 


<*Ip_  (a)]  . 
dz  1 


p(z)g(z) 


In  general,  hydrodynamic  computer  codes  calculate  dp/dz,  the  pressure  gradient, 
by  finite  difference  methods  which  replace  dp/dz  by  Ap/Az.  This  approximation 
results  in  unwanted  accelerations  which  can  be  written  as 


ihz/  \Az/ 
z  z 


For  example,  given  a  code  in  which 


and  assuming  a  scale  height,  H,  constant  over  the  interval*  Az,  so  that  the 
pressure,  p(z),  can  be  written 


p(z)  -  p(zQ)e 


one  can  expand  (Ap/Az)  in  powers  of  Az/H  so  that  to  second  order,  the  net 
accelerations  are 


*<«>  (£) 


,i 


( ,  l: 

•:?v 

-Nil 


For  all  practical  purposes,  this  acceleration  is  negligible  until  the  zone  size, 
A z,  becomes  greater  than  a  tenth  of  a  kilometer.  For  example,  at  an  altitude  of 
250  km,  the  scale  height,  H,  la  about  5  kilometers  and  g  about  900  cm/sec2.  At 
250  km,  the  seals  height  has  its  smallest  value  in  the  atmosphere  so  this 
example  will  give  a  measure  of  the  worst  acceleration.  If  Az  is  0.1  km,  the 
above  relationship  gives  an  acceleration  of  0.05  cm/sec2.  If  Az  is  1  km,  the 
acceleration  is  5.0  cm/sec2. 

The  pressures  and  densities  given  by  the  atmospheric  models  agree  with 
tabulated  values  in  references  1  and  2  to  one  part  in  10*  and  temperatures  to 
ont  part  in  100. 
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APPENDIX 
SUBROUTINE  ATMOS 

The  FORTRAN  subroutine  appearing  In  this  appendix  la  the  annual  wean 
temperate  atmospheric  model.  By  Inputting  an  altitude,  TTY,  the  subroutine 
will  return 

WSP  the  pressure  at  that  altitude 

WSR  the  density  at  that  altitude 

WST  the  temperature  at  that  altitude 

WSI  the  specific  internal  energy  at  that  altitude 

WSU  the  radial  velocity  »  0 

WSV  the  axial  velocity  *  0 

GMONE  y  -  1. 

For  the  annual  mean  tropical  atmospheric  model  or  the  summer  subarctic 
atmospheric  model,  replace  the  data  in  the  TABZ,  TAEL,  TABT,  and  TABP  blocks 
with  the  data  appearing  in  tables  III  or  IV  as  appropriate. 


i 


7c 


SUROOuT  f  Me  A  TWO'i r  Y  .W'lP,  WSP.  WST  »  WS  !  •  W5U •  WS  V «  G*«ONe  ) 
r  i  f'*CCNr  l  Of  J  YA0?t22;  «  TAf?Lf22)  ♦  TART  (22)  •  TABP(?2) 

CALCULATE  ATMOSPHERE 
AxOAD!'!"  Of*  THE  PfiSTH  IN  CM. 

G=ACCFLr  c’Ar  >0N  HUE  TO  GRAVITY  IN  CM . /SFC . /SFC • 
aORC.AftB  FOR  TWO  tFMPEPATF  ATMOSPHERE 

roo  the  tropical  atmosphere 
e'ioj  ,r»t  1  FOR  The  sup  arc  t  J  C  ATMOSPHERE 
Ox CAE  r ONETANT  IN  ifCGS/MOLE/DFG. 
wsmolefl'lap  weight  or  air 


r 

c 


r 

c 


A  =6»  OC-H-7E  +  0P 

r,  =  o,flfwS6e!F+02 

T | fl 4F  +  07 
wsr • ®o6A4F+0 1 

THE  FOLLOWING  DATA  IN  TAB? «  TABL.  T  ART  *  AND  TAPP  APE  FROM  TABLE  2 
FOP  THc  TEMPERATE  ATMOSPHERE.  FOP  THE  TROPICAL  OP  SUBARCTIC  ATMO- 
eoHPOfc  cuaSTITiiTF  ThF  DATA  FPOM  TABLFS  .0  AND  A  PFS.PFCT  I VFLY. 


TAP?  If  THE  BASF  ALTITUDE  OF  EACH  ALTITUDE  GROUP  IN  CM. 
TAB7I  l  J  =0. 

Tio?(  ?  )  a  !  .  1  0  1  nF  +  O-s. 

TA«7C  **  V  =  7.0063E-f06 
TAF/(  4 ) =0.21 62E+06 
T  A  B  ? (  5 ) =4. 7350E  +  06 
T A°7  (  f  j  =*=.?42or+oft 
Tjn7(  7> =6. I 5R1E+06 
TAP? (  °) =7.R«O4E+06 

TAB?!  O)=P,0E+P6 
TAR7( j  0  )  =  1  0  • 0E+06 
T  AP7 (  I  1 )  =  1 1 .0E+06 
T  A°? I  l7)  =  12»0E+06 
TAP?( 1 7) = 1B.0E+06 
TfP?( i 4 ) = 1 6.0E+06 
T  Ap7 ( | P) *  1 7# 0E+06 
ta°7( I  6)  =  1O.0E+06 


TA07(  ,  -7,  =?0.OF+06 
TAPTf ,p) xOO.OF+06 
T£P7(  t^ii^O.OF+Ofi 
TAE7I  -’O)  =GO.0E+O6 
TAp?( 5 1 ) ~60 • OF +06 
TAP7(V>)r7O,0F+06 
f«aL  ic  THF  MOLECULAR  SCALE 
r.r'-'u0  *n  n^r.. /cm. 

TA°Lf  !  >  I  7A7F-OS 

TfiPLf  75=  R.7BC401 77E-0P 
TAt>L(  "»  i  -s  R.  PA2P4P  1  fiF-06 
T  ARL  f  4)=  2.770F0370E-05 
T  Ar’L  (  c  )  =-l  .72246B73E-07 
Tftf»L<  f-  5  =  -  I  •  oepoopopp-05 
ta»L!  ** )  *  --*.n  t  *QT>77ftE“0S 
TAO-Lf  F>=  1  .6CF731 R6F-07 


temdepatupf  goadient  of  each  altitude 


v  \  -n 

xs-v*  -rj  i  -‘ 


;  ; 
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T  AHL.  ‘  *i) 

- 

2  . OR  I  66  '  ML  — OF 

tspl  <  rr> 

0  }  -  2F-05 

T t,  R(_  Mi) 

0.0  7  762  V>'S~0S 

tapl <  i :  ) 

- 

2.0010BB06E-04 

T  API..  {  j  3  ) 

•- 

!  «  4 O RPCi,  .  p . ; T  — 04 

TAPL ( ! 4 ) 

- 

i  .  0040740 ! P  —  04 

TARL C 16) 

-• 

6.0  7F,O6’50  3F— 0  5 

TAPL r ! 6 5 
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